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Light carries angular momentum, and as such it can exert torques on material objects. Applications of these 
opto-mechanical effects were limited initially due to their smallness in magnitude, but later becomes 
powerful and versatile after the invention of laser. Novel and practical approaches for harvesting light for 
particle rotation have since been demonstrated, where the structure is always subjected to a positive optical 
torque along a certain axis if the incident angular momentum has a positive projection on the same axis. We 
report here an interesting phenomenon of "negative optical torque", meaning that incoming photons 
carrying angular momentum rotate an object in the opposite sense. Surprisingly this can be realized quite 
straightforwardly in simple planar structures. Field retardation is a necessary condition and discrete 
rotational symmetry of material object plays an important role. The optimal conditions are explored and 
explained. 

The negative optical torque (NOT) phenomenon manifestly conserves angular momentum 1 4 . Let us consider 
an object illuminated by an electromagnetic wave carrying a z-component angular momentum of L'" c . After 
being scattered by the object, the total z-component angular momentum of the light is V"'. If the con- 
sequence of scattering is such that L s ^ a >L'" c , then angular momentum conservation requires that the difference 
jjm _jjnc must j-, e balanced by a NOT acting on the object. Consequently, negative optical torque is consistent 
with total angular momentum conservation. The question is then whether such phenomenon, which does not 
violate conservation laws, can be realized in practice. If so, is NOT easily observable? We will show that NOT is in 
fact ubiquitous, as it can be realized in many transparent or weakly absorptive structures, in discrete rotationally 
symmetric or even irregular structures, and also in planar or even 3D structures. We shall show that field 
retardation is necessary for achieving NOT and discrete rotational symmetry enhances NOT. The realization 
of NOT adds new possibility to the existing novel 5-12 and practical 13 " 21 technique of optical rotation. It may open 
up new applications in areas where optical micromanipulation play a role 22,23 . We note that although negative 
optical torque here is indeed the angular counterpart of optical pulling force 2427 , it is definitely not a straight 
forward corollary. For example, symmetry is not a key factor in achieving optical pulling force but it proved to be 
of considerable importance in realizing NOT. Also, NOT is ubiquitous while optical pulling force can only be 
observed with a carefully chosen combination of beam profile and particle properties. 



Results 

Explicit numerical demonstration of negative optical torque. Explicit examples of structures where NOT can 
be realized are shown in Fig. 1. For the ease and accuracy of computation and interpretation, we consider planar 
structures composing of dielectric microspheres located on the xv-plane. The axis for which the torque is to be 
calculated is parallel to the z-axis and goes through the center of mass marked in Fig. 1, where some of the 
structures possess discrete rotational symmetry but some do not. We stress that our conclusions are valid not only 
for discrete system comprising of spheres, but also to a continuous piece of material. 

The time-averaged optical torque T acting on the concerned structure is computed by (see Methods) 

(1) 



r=| (^x(T)Wi,is. 

where the time averaged Maxwell stress tensor is given by 2! 



T = 



EE* + c 2 BB*--(|E| 2 + c 2 |B| 2 )l 



(2) 
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Figure 1 | Examples of planar structures that can realize negative optical 
torque, (a)-(g) structures composed of an array of identical dielectric 
spheres. R (a) is the structure (sphere) radius. The crosses mark the centers 
of mass, which coincides with the axis of rotation, (h) shows the coordinate 
system. The incident plane wave has a k- vector (k) along the z-axis. 



the origin of the coordinate system is chosen to be the center of mass 
marked in Fig. 1, a is any surface that encloses the entire structure, 
and r is the position vector originating from the center of mass. 
Irrespective of the structural orientation and morphology, the z- 
component optical torque calculated by Eq. (1) can be written as 



^ hw sca ' 



(3) 



where oj is the angular frequency of incident wave, and m, and m 
characterize, respectively, the azimuthal channels for the incident 
and scattered waves. W^f is the energy extinction rate for incident 
wave from channel m f , and is the rate at which energy is 

scattered into the outgoing channel m. The derivation and the expli- 
cit expressions of and are given in the Methods section. 
As long as the scattering problem can be treated accurately, Eq. (3) is 
exact within the framework of classical electrodynamics. In the 
numerical results shown below, we will consider objects composing 
of an array of spheres as our prototypical system because the scatter- 
ing problem can be solved numerically using the generalized Mie 
theory for multi-sphere to a very high precision. This removes any 
doubt that the NOT effects observed are numerical artifacts. And 
unless stated otherwise, the incident wave is a left circularly polar- 
ized 28 plane wave propagating along the z-axis, with a positive angu- 
lar momentum and an intensity of 1 mW//im 2 . 

Figure 2 shows the calculated optical torques versus ka for three 
morphologically different planar structures with different types of 
incident waves, where k is the wave number and a is the radius of the 
comprising spheres. The results clearly show that NOT can be rea- 
lized in all three structures, which include planar triangular or hexa- 
gonal structures consisting of silicon or polystyrene microspheres, 
illuminated by a circularly polarized plane wave, Gaussian beam, or 
Laguerre Gaussian beam. The optical torques can take negative 
values at some intervals of ka. The three structures in Fig. 2 are not 
of special design. Other planar structures can serve the same purpose. 
Figure 2 illustrates that NOT can be quite universal: it can be 
observed in different structures (planar or non-planar) with different 
forms of incident wave, including plane wave, Gaussian beam, and 
Laguerre Gaussian beam. For the Laguerre Gaussian beam, the smal- 
ler structure (yellow line) do not experience a NOT as the beam 
cannot probe the structural details. While for the larger structure 
(blue line), NOT is still seen. 

Analytical demonstration of negative optical torque and the role 
of retardation. In addition to numerical demonstration, formal 
analytical understanding of NOT can be developed by considering 
the following prototypical system: a pair of non-absorbing dipolar 
spheres located on the xy-plane, illuminated by a z-propagating 
plane wave with arbitrary polarization. Following Ref. 24 and 
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Figure 2 | Optical torque for different structures under different incident 
waves carrying positive angular momentum. Red: Optical torque versus 
varying k for Fig. 1(e) (s r = 2.46, a = 0.49 /im, R= 1.0 \im) illuminated by 
a left circularly polarized plane wave of intensity 1 mW/(im ! . Green: 
Optical torque versus varying a for Fig. 1(a) composed of silicon spheres (s r 
= 12.6096, R = 1.16a, X = 1064 nm) in water illuminated by a left 
circularly polarized plane wave of intensity = 1 mW/(im 2 . Black: Optical 
torque versus varying k for Fig. 1(a) (e, = 2.4649, a = 0.5 /im, R = 
0.58 ;im) illuminated by a left circularly polarized Gaussian beam 
(numerical aperture = 0.9 and total power = 0.05 W). Yellow: Optical 
torque versus fcfor Fig. 1(a) (e r = 2.4649, a = 0.5 \im, R = 0.58 \im) 
illuminated by a linearly polarized 10 mW Laguerre-Gaussian beam with 
topological charge of 1 and a N.A. of 0.9. Blue: Optical torque versus k for 
Fig. 1(e) (s r = 2.4649, a = 0.49 /.im, R = 1.0 /im) illuminated by a linearly 
polarized 10 mW Laguerre-Gaussian beam with topological charge of 1 
and a N.A. of 0.9. Negative torque (T z < 0) means that the light beam 
induces the object to rotate in the opposite sense as the angular momentum 
of the incident photon. 

retaining only the electric dipole term, one can derive the formula 
for optical torque: 



r z = — Re(^| 



where 



/'r 



p±- 



Re LE| E* ± ) Re ( 1 + afi ± + afi« 



+ (\E + \ 



)Im(a/i|i —oifij) 



(4) 



2e lkR 
4m 0 R 3 

e ,kR 

' 4m 0 R 3 



(1-ifcR), 



(k 2 R 2 + ikR-l) 



(5) 



are the eigenvalues of the dyadic Green's function 29 , R is the 
separation between the dipoles, and a is the dipolar polarizability. 
£|l and E ± are the projection of the electric field parallel to or 
perpendicular to the particle axis (i.e. the amplitudes of linear 
polarization), E + and £_ are the amplitudes of left and right 
circularly polarized components. The effect of field retardation 
manifests itself as a factor of e' kR in the expression of wi and /J, ± , 
leading to oscillatory behavior of toque versus kR. We shall see that it 
is this oscillation due to retardation that opens up the possibility of 
NOT. 

The first term in Eq. (4) is proportional to Re(£||£^). It is an 
aligning force, which vanishes for circular polarization. For linearly 
polarization, it vanishes only when either £n = 0 or E± = 0, corres- 
ponding to the case of polarization being perpendicular or parallel to 
the particle axis, respectively. Accordingly, depending on kR, this 
force tries to align the particle chain to a direction that is either 
perpendicular or parallel to the polarization. 

The second term in Eq. (4), which is proportional to (|£+| 2 — 
|£- 1 2 ), also has a simple interpretation: the sign of the corresponding 
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torque is determined by the competition between the amplitudes of 
the left and right circular polarization. Consequently, it is a rotating 
force that will keep the object turning. Depending on kR, a left 
circularly polarized light can induce a positive or negative torque. 
Consequently, this term is responsible for NOT. Moreover, the term 
vanishes identically for linear polarization. 

The two terms in Eq. (4) are competing with each other. When the 
incident wave is predominantly linearly polarized, the first aligning 
term will be dominating and the particle chain will be aligned. On the 
other hand, for circular polarization, the second turning term will 
dominate and the particle chain will be set into rotation. Equation (4) 
contains many salient features of light induced torque. For circular 
polarization, it becomes 



r«=- 



■ (\E+ 1 2 - 1£_ | 2 )/! [Ima+ |a|sin(2fcR+<5)], (6) 



where F x 



and 



8 — tan" 



0 by symmetry, 

k 4 R 4 + 3k 2 R 2 + 9 



32n 2 R 6 4 



(6k 3 R 3 -l8kR)Re{oi} + (9- 
(9-15k 2 R 2 + k 4 R 4 )Re{ct}- 



-l5k 2 R 2 - 
- (6k 3 R 3 - 



-k 4 R 4 )Im{x} 
18kR)lm{a} 



(7) 



(8) 



The first term inside the square bracket of Eq. (6) corresponds to a 
positive bias, whereas the second term carrying sin(2fcR + <5) oscil- 
lates between positive and negative values due to the spatial field 
retardation associated with a finite value of kR. Since |a| > Im a, 
the amplitude of oscillation is larger than that of the bias, accord- 
ingly, the optical torque can sometimes be negative due to spatial 
retardation induced oscillation. This simple exercise not only 
demonstrates analytically the possible existence of NOT, the results 
also show that field retardation is essential to realize NOT. 

For very small separations between the dipoles and a left circularly 
polarized incident wave, Eq. (6) asymptotically approaches 

9Ima 



32rc 2 £ 2 ,tf 6 ' 



(9) 



Since Im a. > 0 by optical theorem 30 , F z > 0. This result is not 
surprising. In the long wavelength limit (/cR<Cl), the wave cannot 
probe the details of the structure. Therefore the sign of the torque 
must be the same as that of the incident wave. In fact, this conclusion 
holds true for any small structures, as we shall see in the following 



Physical origin of the negative optical torque. We now analyze the 
physical origin of NOT for more general structures. Equation (3) 
allows for a heuristic physical interpretation: W^/Hco is the 
extinction rate for incoming photons from channel m v and each of 
these photons carries mfi of angular momentum, therefore the rate 
at which the object acquires angular momentum from all azimuthal 
channels of the incident beam is given by the first term in Eq. (3). 
Similarly WJJJ / hco is the rate of photons being scattered into channel 
m, and each of these photons gives a recoil angular momentum of 
— mh, therefore the rate at which the object obtains angular 
momentum due to recoil is given by the second term in Eq. (3). 
Equation (3) also shows that the transfer of angular momentum 
can be interpreted as occurring in two steps. First, the photons are 
intercepted by the scatterer. This induces an extinction torque given 
by the first term in Eq. (3). For an incident wave consisting of solely 
positive m, components, this extinction torque is positive definite. So 
if NOT exists at all, it must come from the recoil torque in Eq. (3). 
Depending on the material and the structural properties of the object, 
some of the intercepted photons will be absorbed, while some will be 
scattered. The scattered photons exert a recoil toque on the object, as 



given by the second term in Eq. (3), where m < 0 (m > 0) represent a 
positive (negative) contribution. If the negative contribution 
dominates, we will have a NOT. 

The role of symmetry. Symmetry is the explanation to why the 
negative optical torque effect can actually be quite universal. Light 
cannot exert an optical torque on a transparent object possessing 
continuous rotational symmetry, since photons scattered by such 
an object preserve their angular momentum. From this point of 
view, the phenomenon of optical torque is a consequence of 
broken rotational symmetry. NOT is therefore inevitably tied to 
the symmetry of the structure and incident wave. 

To see this, it suffices to consider objects with discrete rotational 
symmetry of degree m s and for an incident wave with a single azi- 
muthal channel m, > 0. In wave scattering theory, based on irre- 
ducible tensor approach, wave can be decomposed into modes 
characterized by two mode indices: total angular momentum index 
Z, termed also angular momentum channel index, and z-component 
of angular momentum index m, termed also azimuthal channel index 
in this paper. We will see that it is discrete rotational symmetry that 
allows NOT to be observed ubiquitously. 

The angular momentum selection rule for photons after being 
scattered by such object is 

(10) 



m — W,- + n x m.. 



where m is an allowed channel for the scattered wave and n = 0, ±1, 
±2 and so on 31,32 (see part 1 of Methods for a proof). Equation (10) is 
the angular counterpart of the grating effect for a periodic surface 
illuminated by a plane wave. Although Eq. (10) is analogous to linear 
grating (periodic in space), the rules are correct for all structures in 
the rotational case because any structure is invariant under a rotation 
of 2tl With Eq. (10), the z-component optical torque of Eq. (3) can be 
re -written as 



mfi 

— W , 
not 



(mi) 



^(m, + nxm s )-H*' + " x "'». (11) 



In principle, the positive and negative n channels contribute to NOT 
and positive optical torque (POT), respectively, while the n = 0 
channel, representing no change in the angular momentum of the 
photons, produces no torque. However, the contribution for chan- 
nels with \m\ 2> kR can be rather weak. When \m\ 2> kR, the cen- 
trifugal effect forces the amplitude of the mode m to decrease rapidly 
for region where R < |m|/fc 33 . There is very little overlap between the 
mode m and the object, forbidding the object to scatter wave into 
these channels. Consequently, these channels contribute negligibly 
to optical torque. 

If m s is very large, the condition of |m| 2> kR will be met for all 
channel except the n = 0 channel, even for the low order modes with 
n = ±1. In fact, continuous rotational symmetry is an idealization 
rather than reality. An object with very fine sub-wavelength details 
(such as a small object with a large m s ) that is beyond the diffraction 
limited resolving power of the incoming wave can be thought to have 
continuous rotational symmetry as far as light scattering is con- 
cerned. In this case, the contributions from all channels become 
vanishingly small except the n = 0 channel. Therefore the wave 
exerts vanishingly small torque on such object with high discrete 
rotational symmetry and relatively small size, just like a structure 
with continuous rotational symmetry. The incident wave m, can only 
be scattered into the azimuthal channel characterized by the same m,. 
Equation (3) thus reduces to T z = mfiW^^/hco, where 



w ( , m,) = wi , 

fv abs " v ext 



W^' is proportional to the energy absorption rate 



in channel m,. Since both W%£~ and m, are positive definite, it follows 
that r„ > 0. 
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In the other extreme situation of m s <?Cm iy according to Eq. (3), the 
torque is positive definite as 

\ " / (12) 

(W ext -W sca ) m { Wai, s 
= rtii = > 0, 

ft) CO 

where W a ^ s = (W ext — W sca ) is the rates at which energy is absorbed. 

If NOT can exist, it must be for some moderate value of m s < kR so 
that the corresponding partial waves can reach the scatterer to trans- 
fer angular momentum. This is in fact a manifestation of the optical 
diffraction limit. If the small structure has some very fine sub- 
wavelength details, such as patterns with high degree of rotational 
symmetry located in a small region, the incident wave will not be able 
to probe them. Thus the information on the symmetry will be lost, 
eliminating the possibility of observing NOT. In general, NOT can be 
easily observed if m s < Z max , where l max = kR + 4.05(fcR) 1/3 + 2 is an 
empirical formula for the angular momentum truncation order of 
the scatterer 31,32 . 

Next, we argue that a larger m s favors NOT as long as \m\ is not 
significantly larger than kR. In Fig. 3, the optical torque acting on 
objects composing of the same type of spheres but with different 
values of m s (m s = 1, 2, 3) are shown. By comparing m s = 1 to m s 
= 3, one sees that NOT is more easily observed for a larger value of 
m s . These can be explained as follows. We note that according to our 
simulations, and also Eq. (6), the total torque consists of a positive 
bias part and an oscillatory part. It is the oscillatory part that gives 
NOT. For the object with m s = 1 (red line), the torque acting on all 
three spheres are different. Therefore the negative torque acting on 
one sphere may be canceled by the positive torque acting on another. 
This weakens the oscillatory part of the total torque and reduces the 
possibility of NOT. For the object with m s = 3 (blue line), the torque 
acting on each sphere is the same by symmetry, so the torques rein- 
force each other rather than cancel each other. This preserves the 
strength of the oscillatory part of the total torque and allows NOT to 
be observed more easily. More discussion about the role of symmetry 
can be found in the Supplementary Information. 

Discussion 

For the more complex systems we have considered (arrays of spheres 
as shown in Fig. 1), we further argue here that spatial retardation is 
also essential for NOT, as in the simple example we discussed in the 
previous section. We show in Fig. 4 the optical torque due to indi- 




Figure 3 | Optical torque versus ka for objects with different discrete 
rotational symmetry. The objects are composed of identical dielectric 
spheres (e t = 2.4649). The red, black, and blue lines correspond to m s = 1, 
2, and 3, respectively. The degree of discrete rotational symmetry is seen to 
significantly enhance negative optical torque (i.e. higher m s favors negative 
torque). 



vidual azimuthal channel for structures with three-fold (m s = 3) or 
six-fold {m s = 6) rotational symmetry. According to Eq. (10), for 
three-fold symmetric structures showed in Fig. 4(a)-(b), light is scat- 
tered into azimuthal channels characterized by m = In + 1, where n 
= 0, ±1, ±2, ... etc. The m = 1 channel [i.e. n = 0 in Eq. (10)] does 
not contribute to torque for lossless spheres. Light scattered into 
azimuthal channels m = —2, —5, — 8 and so on will contribute to 
positive optical torque, whereas the channels characterized by m = 4, 
7, 10 and so on will induce NOT. As multiple scattering is not strong 
in dielectrics, the recoil torque is dominated by n = ±1, correspond- 
ing to m = ntj ± m s = — 2 or 4. In Fig. 4 the oscillatory behavior of the 
torques due to individual partial waves versus wavelength and struc- 
tural size is the signature of retardation, while the peaks for n = ± 1 
channels correspond to Mie type resonance. It is precisely due to this 
retardation induced oscillation that the total torque takes positive or 
negative value. The situation for the six-fold symmetric structures 
showed in Fig. 4(c)-(d) is similar. 

Finally, we argue that small object does not favor NOT. Consider, 
for example, a small three-fold {m s = 3) rotationally symmetric 
lossless structure illuminated by a left circularly polarized plane wave 
(m,- = 1). If the structure is small such that wave components with 
higher angular momentum cannot reach the scatterer and the trun- 
cation angular momentum index Z max = 1, then the torque is 
approximately equal to zero since only n = 0 term appears in Eq. 
(1 1). As the size of the structure increases such that we have to take 
'max = 2, the second term in the series (corresponding to n = — 1) is 
now significant, which is given by (see Methods for the definition of 
partial torque z m that is due to each individual outgoing azimuthal 
channel m) 

T- 2 = m s W_2,2/(»- (13) 

Equation (13) is positive definite and therefore corresponds to POT. 
One then expects that a small scatterer cannot experience a NOT. If 
we consider larger particle, NOT will become possible. For example, 
the third term in the series (corresponding to n = 1) is given by 

t 4 = — m s Wi^/co, (14) 

which, when larger in magnitude than x_ 2 in Eq. (13), results in 
negative total torque. 

In summary, we showed that negative optical torque is ubiquitous 
in simple configurations and can be induced by a broad class of 
incident waves. We note that while negative optical torque can be 
viewed as the angular counterpart of the optical pulling force, the 
conditions under which the "negative response" can be realized are 
very different. The optical pulling force is achievable only for special 
kinds of optical beams acting on some specific particles 24 27 whereas 
NOT can be realized rather easily. NOT can occur if the recoil torque 
happens to be negative and larger in magnitude than the extinction 
torque, which is made possible by retardation effects and enhanced 
by discrete rotational symmetry. We remark that our theory is valid 
for both spin and orbital angular momentum (i.e. arbitrary m,) of 
light. 

Methods 

Derivation of the azimuthal mode selection rule in electromagnetic scattering. We 

shall prove that when an incident wave characterized by a single azimuthal number 
m; 30,32 impinges on a structure with m s -fold rotational symmetry, it will be scattered 
into azimuthal channels given by the azimuthal mode selection rule in Eq. (10). 

Let us start with the total electric field for a general scattering problem, which is 
given by the Lippmann-Sch winger integral equation 34,35 

E(r) = E inc (r) + |G(r,r')- [ £ ,-(r')- l]/: 2 E(r')dr', (15) 

where E (J1L (r) denotes the incident electric field, £ r (r') is the position dependent 
relative dielectric constant, and k = w/c is the wave number in the background. The 
second term on the right hand side of Eq. (15) gives the scattered field: 
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Figure 4 | The contribution (in absolute value) to the torque by individual partial wave azimuthal channels versus kR. (a)-(b) three-fold and (c)-(d) 
six-fold rotationally symmetric dielectric structures (e r — 2.4649, a — 0.49 firri). The incident field is a left circularly polarized plane wave (intensity is 
1 mW//im 2 ). R — 0.58 (im for (a) and R — 1 urn for (b)-(d). The oscillations are due to spatial retardation. 



E sca = J G(r,r')-[ £r (r')-l]rC 2 E(r')dr', (16) 
where G(r,r') is the free-space dyadic Green's function 29 , which can be written as 
G(r,r')= *«rmO« m i(rM)0^m)i(r',& , 4'), (17) 

<x,l,m 

where O amn represents vector spherical wave functions M mn> N mn and L m „ 29,36 for a = 
1, 2, and 3, respectively. The ^-dependent part of O am j (r, 0, <j>) is of the form e' m ^ and 
similarly O a (- m )/ (r', 9', <f>') depends on <j)' through e~ im,i '' . As a result, each term 

characterized by m in the dyadic Green function G(r,r') given by Eq. (17) shows 
e im(<i> -•!■') azmiuma i dependence. 

Inserting Eq. (15) into Eq. (16), we arrive at 
E sca - jG(r,r')-[^(r')-l]fc 2 E inc (r')dr' 

+ |G(ry)-G(^01^W-l][^(0-l]fc 4 Einc(rWdy (18) 

where ■ ■ ■ denotes higher order scattering terms. Note that the azimuthal dependence 
of the incident wave E iMC (r') on <f>' is given by e' m '^ , while for a structure with m s -fold 
symmetry, the dielectric constant may be expanded in a Fourier series 

CO 

e r (r')-l= YJ A » e ™''' ( 19 ) 

n= — oc 

Therefore, for the first term on the right hand side of Eq. (18), the integration over 
azimuthal coordinate (j) ' reads 

/■2tt 

YJ • • • e - im *' ( VJ Aj nm -vy m >* df , 



which requires m — m ( + nm s in order to have a non-zero contribution to E sca . 
Similarly, for the second term on the right hand side of Eq. (18), the integration over 
(j)' and <f>" can be written as 

It then follows straightforwardly that the non-vanishing terms should satisfy m ' = 
rrij + n'm s and m = m' + nm s> implying that only azimuthal modes with m = (n + 
n')m s + m, survive in the scattered field. 

The procedure applies straightforwardly to higher order scattering terms in Eq. 
(18), leading to the conclusion that the scattered field contains only azimuthal modes 
with m satisfying m = m, ■, + nm s , where n is an integer. This is actually an angular 
analogue of the much better known phase matching condition for a surface of peri- 
odic lattice derived from Floquet's theorem. 

The contribution of individual azimuthal channel to optical torque. The time- 
averaged optical torque (T) t acting on a scatterer is obtained from the surface 

integration of the time-averaged angular momentum flux tensor ^K^ by 

(r) f =-J n(K) t dS=-JJ n-((T) ( xr)dS, (20) 

where a can be any surface that encloses the entire scatterer, and n denotes the unit 
outward normal at surface o. The time-averaged Maxwell stress tensor ^T^) is given 
by 

(f>, = iRe[ f: „E e E; + , lo H e H: - ^(F. 0 E e -E' e + t i 0 H e -H' e )i}, (21) 

where e 0 and /i 0 are, respectively, the permittivity and permeability in free space, and I 
is the unit tensor. The electric (magnetic) field E c (HJ is a sum of incident field E, (H ; ) 
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and scattered field E s (H s ), namely, E e = E, + E s and H e = H, + H,. These can be 
expanded in terms of vector spherical wave functions (VSWFs) 37,38 : 

00 I 
1=1 m=-l 



C0 '' 1-1 m=-l 
GO / 



(22) 



where Z max is the truncation order for the angular momentum index in the infinite 
series. Equation(26) implies that all partial torques associated with positive (negative) 
n contribute to NOT (POT). The partial torques z m for some different allowed m s 
(with n = ±1, ±2, ±3) are shown in Fig. 4 for the case where a left circularly plane 
wave {with m, — 1) impinges on scatterers with m s — 3 and 6. From Fig. 4, it is found 
that the absolute values of partial torque due to the same |n| are quite close, since they 
contain roughly the same number of terms according to Eq.{26). In particular, all 
terms show oscillatory behavior versus kR due to field retardation. It is this retarda- 
tion-induced oscillatory behavior that gives rise to the competition between partial 
torque with positive and negative n, resulting in NOT in some intervals of values in 
kR. 



H s = — Yl M^N2(M) + fl*iM2(*yr)]» 

0) /* 1=1 ,„=-! 



where E m 



■ and y n 



2/+1 (l-m)\ 



The expansion coefficients p mi and 2 



l{l+l)(l+m)\ 

1m\ ( a mi an< i &ml) are known as partial wave expansion coefficients for incident 
(scattered) waves. Here, p m \ and q m i can be obtained from an overlap integral between 
the incident field and VSWFs. For a collection of spheres a m \ and b m \ can be obtained 
by solving the scattering problem using the generalized Mie theory for multispheres. 

By the orthogonality of the VSWFs, the integration in Eq.(20) can be worked out 
analytically {for an arbitrary incident field) at a surface far away from the scatterer to 
give the z-component T z of the time-averaged optical torque, 



2kf,\E[ 



fc 3 



;=i m =-\ 



On the other hand, with the incident and scattered field expanded in terms of VSWFs, 
the rate W sca at which incident wave energy is scattered and the rate W ext at which the 
energy is intercepted can be written as, see, e.g., Ref. 38 



where 



m 

W m =-£n-(P„,),dS=E^ ) > 



WJS- 7-, l^(\ami\ +\b m {\ ) 



(24) 



k- 



2?CSocLEo 



k 2 



and the timed-averaged Poynting vectors {P sm ) t and (P exf } ( are given by 



(P ra ) f = -Re[E s xH*], 



(p e x.>,=2 Re [ E i xH r+^ xH ri' 



For an incident wave from a single azimuthal channel m„ the partial wave 
expansion coefficients satisfy p m \ =p m \& mm . and q m i — q m iS mmt , where 3 mm > is the 
Kronecker delta. So Eqs. (23) and (24), together with W ext = W sca for lossless scat- 
terer, imply that 



r z = - 



lne a LEn 



E mi (l 0 >" 



E' 



27t£o|£o|' 



^2(mi-m)(\a ml \ 2 + \b,„ 



where W m \ — 



27t£ 0 c|£ 0 



=-~ E » E w >»<< 

«=— CO I = max{|m|,l} 

■ (|«m;| 2 + |^m/| 2 ) >0 with Wsca = W m |, and m - 



(25a) 



(25b) 



(25c) 



nm s for a scatterer with m s -fold rotational symmetry, see, Eq.(lO). The first sum- 
mation in Eq. (25a) corresponds to the extinction torque, which has the same sign as 
m h leading definitely to a POT, while the second summation in (25a) corresponds to 
recoil toque that can be positive or negative, since m( = m, + nm s ) and m\ may be of 
the same or opposite signs. Equation (25c) indicates that in order to exert an optical 
torque on a lossless object, the latter must scatter the incident wave into azimuthal 
channel m that differs from the incoming one m, (such that n 9^ 0). Based on Eq. (25c), 
one can separate the optical torque F z into partial torque z rn that is due to each 
individual azimuthal channel m of the scattered field 



E-=E^ 



l = mi\x{\mi + nm s \,l} 



(26) 



— ^ ^2 mRelp^a^^q^b^-ija^ 2 ^\b m! \ 2 )}. (23) 7. 
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